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ABSTRACT 

We apply a spherical harmonic analysis to the Point Source Redshift Survey (PSCz), 
to compute the real-space galaxy power spectrum and the degree of redshift distortion 
caused by peculiar velocities. We employ new parameter eigenvector and hierarchical 
data compression techniques, allowing a much larger number of harmonic modes to 
be included, and correspondingly smaller error bars. Using 4644 harmonic modes, 
compressed to 2278, we find that the IRAS redshift-space distortion parameter is 
P = 0.39±0.12 and the amplitude of galaxy clustering on a scale of k = 0.1 /iMpc -1 is 
Agai(O.l) = 0.42 ±0.02. Combining these we find the amplitude of mass perturbations 
is A m (0.1) = (0.16 ± 0.04)$1~ 06 . A preliminary model fitting analysis combining the 
PSCz amplitudes with the CMB and abundance of clusters yields the cosmological 
matter density parameter fl m = 0.16±0.03, the amplitude of primordial perturbations 
Q = (8.4 ± 3.8) x 10" 5 , and the IRAS bias parameter b = 0.84 ± 0.28. 
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1 INTRODUCTION 

The extraction of cosmological parameters from surveys has 
entered a new phase, with the advent of very large data sets. 
But the prospect of accurately determining a wide range of 
parameters is offset by the difficult task of manipulating 
these large data sets without loosing important parameter 
information. In the case of analysing near-all sky redshift 
surveys we have developed a method based on the harmonic 
decomposition of the survey into spherical harmonics and 
radial spherical Bessel functions (Heavens & Taylor 1995, 
Ballinger, Heavens & Taylor 1995, Tadros et al 1999; here- 
after HT, BHT and T99, respectively). The use of a spheri- 
cal harmonic decomposition allows the accurate analysis of 
the radial redshift-space distortion effect, without using the 
small-angle or distant observer approximations, and the nat- 
ural inclusion of angular and radial window functions. The 
parameters of interest are the degree of redshift-space dis- 
tortion, parameterised byQ 



(i) 



* Since differently selected galaxies cluster differently, each se- 
lection may have its own bias parameter and /3. In this paper we 
shall use /3 and b to refer exclusively to IRAS selected galaxies. 



where /(fi m ) = d\n8/dlna « f^ 6 (Peebles 1980) is the 
growth rate of perturbations and b is the linear bias param- 
eter defined by 



>(*), 



(2) 



and an estimate of the undistorted galaxy power spectrum, 
Pg a i(fc). This approach combines the spherical harmonic 
decomposition with a mode-by-mode maximum likelihood 
analysis, and has been applied to the 1.2Jy survey, a 1:6 
subsample of the IRAS galaxy survey (HT, BHT) and the 
full 1:1 sample, the PSCz (T99). T99 used the most sophis- 
ticated analysis to date, and found a distortion parameter 
of p — 0.58 ± 0.26 (marginal error) and an amplitude for 
the real space galaxy power spectrum at k — O.lftMpc -1 
of A .i = 0.42 ± 0.03, where A 2 (fc) = k 3 P(k)/2n 2 . In this 
paper we employ new data compression methods to improve 
on this analysis and obtain more accurate parameter deter- 
minations. We describe the methods more fully in Ballinger 
et al (2000). 

The limiting factors in our previous analyses were CPU 
time and stability. To complete the likelihood analysis the 
data covariance matrix of the full data set must be inverted 
at each point in parameter space. Small numerical errors can 
make this procedure unstable; data compression (Section 2) 
is a great help here, as the high signal-to-noise modes are 
well-behaved. The CPU factor can also be an issue when one 
wishes to investigate systematic effects in data sets. This 
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is as much an issue in analysis of the Cosmic Microwave 
Background (CMB) as it is in galaxy redshift surveys, such 
as the PSCz, the 2-degree Field (2dF) or the Sloan Digital 
Sky Survey (SDSS). 

The problem of analysing large data sets was addressed 
by Tegmark, Taylor & Heavens (1997; TTH) who considered 
the question of what was the optimal transformation of the 
data for estimating a given parameter, where the model data 
covariance matrix could be an arbitrary, nonlinear function 
of the desired parameter. The optimal transformation should 
have the properties of maximising the information content 
about the parameter in the minimum number of eigenmodes. 
By ordering modes the ones with the most information could 
be selected and the rest discarded. That the data can be 
ordered this way can be understood if one considers that 
the data may contain noisy or strongly correlated modes 
that add little information about the parameter of interest. 
Choosing the transformed data to have diagonal covariance 
also decreases the computation time. While the data covari- 
ance is only diagonal at one point in parameter space, the 
removal of correlated and noisy data by trimming produces 
a numerically stable inversion of the covariance matrix. 

This procedure sounds similar to Principle Component 
Analysis (PCA) or signal-to-noise eigenvalue analysis (SNA; 
Bond 1995, Vogeley & Szalay 1995), but has important dif- 
ferences. We have previously referred to the procedure as 
Karhunen-Loeve methods, but it is in fact more general, 
so we shall refer to our method henceforth as Optimal- 
Mode Analysis (OMA). In this paper, we introduce two 
new methods for accurate parameter estimation, and refer 
to the whole method as Generalised Optimal Mode Analysis, 
or COMA. 

The additional features of COMA can be split into two 
parts, dealing with multiparameter estimation and the sta- 
ble handling of data compression. In multiparameter esti- 
mation it is useful to note that for highly-correlated mul- 
tiparameter estimates (highly-elongated likelihood surfaces, 
not aligned with parameter axes), the marginal error in the 
parameters is determined by the length of the longest likeli- 
hood principal axis. We therefore want to optimise to keep 
this length as short as possible. This process is called param- 
eter eigenmode analysis, and was introduced in Ballinger 
(1997), with some results being presented in Taylor et al. 
(1997). The second, optional part of COMA is to split the 
original data into subsets, optimising each subset, and then 
combining the best modes together. This procedure can be 
used hierarchically, to reduce a very large number of modes 
to a manageable size. This process we refer to as hierarchi- 
cal data compression. These methods will be detailed in a 
companion paper (Ballinger et al. 2000). 

In this paper we combine our spherical harmonic de- 
composition, parameter eigenmode analysis and hierarchical 
data compression methods to analyse the PSCz. We study 
both nonlinear multiparameter estimation, the redshift- 
space distortion parameter and the amplitude of power, us- 
ing hierarchical data compression, from the harmonic modes 
of the PSCz survey. The increase in analysing power using 
these methods allows us to increase the number of modes 
available for study, and hence a corresponding increase in 
accuracy of our results. 

Padmanabhan, Tegmark and Hamilton (1999) have also 
used the spherical harmonic decomposition to analyse the 



CFA/SSRS UZC galaxy redshift survey, while Hamilton, 
Tegmark & Padmanabhan (2000) have applied the analy- 
sis to the PSCz redshift survey. In both cases they employ 
Karhunen-Loeve methods to estimate the quadratic band- 
power in these surveys. In the former survey the band-power 
was measured in redshift-space, while in the latter analysis 
of the PSCz they measured the real-space galaxy-galaxy, 
galaxy-velocity and velocity-velocity power spectra, and es- 
timated (3 from a least squares fit to the ratios of the galaxy- 
velocity to galaxy-galaxy and velocity-velocity to galaxy- 
galaxy power spectra. 

The paper can be summarised as follows. In Section 
2 we briefly review the spherical harmonic decomposition, 
OMA, parameter eigenmodes and hierarchical data com- 
pression methods. In Section 3 we perform an maximum like- 
lihood analysis of the PSCz, for the redshift-space distortion 
parameter and real-space clustering amplitude of galaxies. 
Our results are presented in Section 4, and conclusions are 
made in Section 5. We begin by a brief review of our data 
analysis methods. 



2 DATA ANALYSIS 

2.1 Spherical Harmonic decomposition 

Following HT, BHT, and T99 we can decompose the density 
field of galaxies in a redshift survey into harmonic modes 

plmn = Clin J d :i s p(s)w(s)j e (ks) Y e * m (fl) , (3) 

where Y{ m , are spherical harmonics, ji are a discrete set of 
spherical Bessel functions, w(s) is an adjustable weighting 
function and s is the redshift-space position variable. The 
cii n are normalization constants and k — fe;„ are discrete 
wavenumbers. Each mode was weighted with a Feldmann- 
Kaiser-Peacock (1994) weight, w(k,s) = [1 + <j>(8)P{k)]~ 1 , 
where cf>(s) is the (redshift-dependent) average number den- 
sity of galaxies in the survey. 

The observed coefficients, Di mn , can be related to the 
underlying linear density modes, 8( mn , by (HT, BHT, T99) 

D = p-p = SW($ + l3V)8. (4) 

The transition matrices S, W , <I> and V correspond to the 
effects of small-scale random radial velocity distortions, the 
angular window function, the radial galaxy selection func- 
tion and linear redshift space distortions, respectively. The 
mean field, po, is nonzero due to the radial selection function 
and angular window function. Note that there is no matrix 
multiplication implied between the angular matrix, W , and 
the radial matrix, <l? + f3V , since these are orthogonal. 

The small-scale radial velocities can be accurately 
modelled by multiplying the galaxy power spectrum by a 
Lorentzian function. This implies that each mode is mul- 
tiplied by the square-root of a Lorentzian. Inverse Fourier 
transforming we find that the density field should be con- 
volved with the function 

S( X ) = 2 -^K (f X ) , (5) 

where K n is an n th -order modified Bessel function and a v is 
the 1-d velocity dispersion. (Note that in the analysis of T99, 
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the 3-d velocity dispersion was used incorrectly. Changing to 
the correct value had little effect. However, as we are push- 
ing our model to higher wavenumbers here, it is more im- 
portant to have the correct value of the velocity dispersion.) 
The transition matrix, S, is found by a spherical harmonic 
transformation of S(x — y). 

Two immediate advantages of this treatment are ac- 
counting for the effects of the monopole and dipole modes. 
The monopole mode contains information about the mean 
density of the survey. In previous methods this can bias 
down the estimated power as the mean is estimated from 
the survey itself and may not be the true mean (Tadros 
& Efstathiou 1996). In our treatment the monopole mode 
can be removed, effectively removing this bias (T99). As the 
PSCz is not all sky, some aliasing takes place at the few per- 
cent level, and we include the effects of this leakage, through 
W . The dipole includes contributions from our own motion 
in the redshift space distortion. Again we can mostly remove 
this by ignoring the dipole mode, and accounting for aliasing 
from the angular mask (T99). 

The distribution of linear harmonic modes can be de- 
scribed by a multivariate likelihood, 

-2hxC{D\j3,P(k)] = D f C X D + lndetC, (6) 

where C = {DD } is the data covariance matrix. Details 
for dealing with non-axisymmetric angular window functions 
in the covariance matrix, as well as further details of the 
likelihood analysis, are given in T99. 

2.2 Generalised Optimal Mode Analysis (GOMA) 

The technical advance which allows us to reduce the error 
bars in the determination of /3 and the power spectrum is 
a new optimised form of data compression, which we call 
generalised optimal mode analysis. Details of the method 
will appear elsewhere (Ballinger et al. 2000), but we sketch 
the main ingredients here. 

The need for data compression is twofold: first, the 
speed of analysis generally scales as N 3 , where N is the 
number of modes considered. These modes might be, for ex- 
ample, spherical transform coefficients. Secondly, since the 
covariance matrix is computed numerically, small numeri- 
cal errors can lead to a non-positive definite matrix, which 
makes no physical sense; even a single negative eigenvalue 
out of several thousand is fatal. Instead of working with the 
full set of spherical transform coefficients up to some maxi- 
mum wavenumber, we form orthogonal linear combinations, 
and use these as the data. GOMA consists of two parts; one 
is parameter eigenmode analysis, where instead of choosing 
/3 and Ao.i as the two parameters to be determined, we 
introduce a new parameter, which runs along the longest 
likelihood ridge of Figure 1, and use the data compression 
methods of TTH to make the likelihood as narrow as pos- 
sible in this direction. Since the length of this controls the 
marginal errors in both j3 and Ao.i, the method is very ef- 
fective. It is worth noting that only OMA will determine 
the best set of eigenmodes following this procedure, since 
we now have a linear combinations of parameters which are 
nonlinear in the covariance matrix. 

The second (optional) ingredient is hierarchical data 
compression. We cannot find the optimal linear combina- 
tions of the entire mode set, because of the numerical prob- 



lems identified above. Instead, we form optimal orthogo- 
nal linear combinations of subsets of the modes (in discrete 
wavenumber ranges), and then combine the best modes to 
form a new set. This process can be repeated hierarchically 
to produce a near-optimal set of modes. Note that when 
mode sets are combined they are not orthogonal and we use 
the correct correlation properties. 



3 APPLICATION TO THE PSCZ 

We have applied our methods to the new IRAS PSCz red- 
shift survey (Saunders et al 1999), with a flux cut of 0.75Jy, 
and the ultra-conservative mask created by T99, correspond- 
ing to optical extinctions of Ab = 0.75mag and excluding 
35% of the sky. Our flux cut is above the formal limit of 
of 0.6Jy for the catalogue, but motivated by a systematic 
effect we found linked to flux in T99. This appeared as a 
flux-dependency on the amplitude of real space perturba- 
tions, where the low-flux galaxies had roughly a factor of two 
higher clustering amplitude, the cause of which we were un- 
able to identify. However by cutting the catalogue at 0.75Jy 
the effect could be removed. This cut in flux, along with the 
more conservative mask left us with around 7000 galaxies. 

The survey was put in a sphere of radius r max = 
200/i~ 1 Mpc and transformed into spherical harmonic 
modes. We analysed all the modes down to k = 0.2/iMpc -1 , 
with n = 1 - 14 and I = 2 - 36 yielding a total of 4644 
harmonic modes. Modes in the range n — — 20 and 
I — — 50 were used for the convolution matrices. Since we 
use a higher cut in wavenumber, k, than in the analysis of 
T99 we must be cautious about introducing nonlinear modes 
which go beyond our analysis. The main concern is the effect 
of "fingers-of-god" contaminating our analysis. The effect of 
these would be to lower the measured value of (3, since its 
effect is to elongate structure along the line of sight, and 
lower the clustering amplitude. We have tested our methods 
using CDM simulations under similar conditions, and find 
that our correction for radial small-scale velocities, equation 
(^), is accurate (Ballinger et al 2000). 

The 4664 modes were compressed down to 2278, using 
the hierarchical compression method once, and a new like- 
lihood analysis applied. In T99 around 1300 modes where 
analysed. In this analysis we have partly made use of the 
increase in computing power in the time between the two 
analyses. However the real limiting feature of our previous 
analysis was numerical instability problems in the covariance 
matrix. Small errors can produce a covariance matrix which 
is not positive-definite, which makes no physical sense (a 
probability distribution in data space which grows exponen- 
tially along one principal direction). A great advantage of 
the compression mechanism is that this numerical problem 
is completely solved: the best modes for parameter estima- 
tion are well-behaved, and the covariance matrix inversion 
proceeds smoothly. Hence the overall time for computing 
our results (around 1 weeks CPU time) remained constant. 



4 RESULTS 
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Figure 1. Likelihood contours for the parameter space of redshift 
space distortions, /3, and real-space galaxy power, Ao.i, measured 
at k = 0.1/iMpc -1 for the PSCz survey cut at a flux limit of 
0.75Jy. The larger set of contours is the results of the Tadros et 
al (1999) analysis using 1300 modes and fc ma x = 0.13 /iMpc -1 . 
The smaller set of contours is the present analysis using 4644 
modes compressed into 2278 modes and fc max = 0.20 Mlpc" 1 . In 
both analyses a conservative mask was used leaving a total sample 
of 7042 galaxies. The contours are spaced by Aln£ = —0.5. 



4.1 Likelihood Analysis 

Figure |l| shows the results of our likelihood analysis for the 
f3 — Ao.i plane. The larger contours are the results of T99 
for 1300 modes. Contours are plotted at intervals of —0.5 
in ln£ from the maximum. The smaller set of contours are 
for our new analysis for 4644 harmonic modes compressed 
to 2278. We have used a CDM-type power spectrum with 
shape parameter T = 0.2 to calculate the covariance ma- 
trices, leaving the amplitude a free parameter, and in the 
wavenumber-dependent mode weighting function. We used 
a value of a v = 224 kms" 1 for the 1-dimensional velocity 
dispersion of galaxies in the scattering matrices, but have 
also experimented with 112 kms -1 and 336 kms -1 . 

The increase in parameter information yields a new, 
lower value of 



(3 = 0.39 ±0.12, 



(7) 



where we quote the marginalized errors. This is consistent 
with our previous result, but with a much smaller error el- 
lipse, around a factor of 3 improvement, which significantly 
rules out both j3 = 1 and /3 = 0. The amplitude of the real 
space galaxy power spectrum is 



Ao.i = 0.42 ± 0.02 



(8) 



which is again consistent with the results of T99 with a 
slightly smaller error. The covariance parameter of f3 and 
Ao.i is r — 0.82, estimated from the ellipticity of the er- 
ror ellipse. We find no evidence that our analysis is being 
biased by nonlinear effects, since then we would expect a 
sudden change in both the value of j3 and a drop in the am- 
plitude of real-space perturbations, neither of which is seen. 
We also find that the maximum likelihood values are almost 
unchanged, moving by less than A/3 « 0.1, if we change 




<=> 0.01 



Figure 2. The real-space power spectrum of IRAS galaxies. On 
large, linear scales the analysis of Tadros et al (1999) is used 
to estimate the real-space band-power in five passbands. On 
smaller scale we have plotted the real-space power spectrum esti- 
mated from the analysis of Saunders et al (1997) from the cross- 
correlation of the QDOT and QIGC surveys. The lighter points 
are the mass amplitudes estimated from the redshift space distor- 
tion parameter. The last point is the mass amplitude estimated 
from the abundance of clusters. The model fit is a ACDM model 
with f2 A = 0.84 and U m = 0.16, h = 0.65. The lighter line is the 
linear fit, while the solid line uses the nonlinear transformation 
of Peacock & Dodds (1996). 



the velocity dispersion to 112 or 336 kms -1 , suggesting that 
nonlinear effects are not significant for this analysis. 

As the amplitude of galaxy perturbations is propor- 
tional to the bias factor, we can combine (3 and the linear 
galaxy power spectrum to estimate the amplitude of the 
mass power spectrum, 



(k) = (3A gal (k)n- 



(9) 



The real-space power spectrum of the PSCz was estimated 
by T99, and is plotted in Figure 2. Also shown is the am- 
plitude of mass perturbations for 52 m = 0.16. The point 
at k = 0.23 /iMpc" 1 is the mass amplitude on a scale of 
8 h~ 1 Mpc, estimated from the abundance of clusters (Henry 
& Arnaud 1991, White, Efstathiou & Frenk 1993, Viana & 
Liddle 1996, Eke, Cole & Frenk 1996), 

-0.47 



as = 0.56S1„ 



(10) 



with a conservative error of around 30%. The mass ampli- 
tude from the PSCz at k = 0.1 /iMpc" 1 is 



A m (o.i) = (o.i6±o.o4)n; 



(11) 



In Figure 2 we plot the mass amplitudes for the PSCz 
and the cluster abundance. There is a clear consistency be- 
tween these two estimates of the mass amplitudes. Assuming 
bias is linear on large scales, estimating mass amplitudes 
from large-scale redshift surveys is simpler, and hence in 
principle more robust, than the abundance of clusters ar- 
gument. In addition this can be used to sample the linear 
spectrum of mass perturbations on a range of scales, rather 
than being restricted to one scale. 



© 0000 RAS, MNRAS 000, 000-000 



Redshift- space distortions of the PSCz galaxy catalogue 5 



4.2 Cosmological Parameters 

Although the shape of the real-space, linear galaxy power 
spectrum can be used to constrain models of structure for- 
mation, the range of points in the linear regime is limited. 
A better test of models is to compare the mass amplitudes 
from the PSCz against the amplitude estimated from the 
CMB. To span the scales between different measurements 
we need to assume a model linear mass power spectrum. We 
shall assume that a standard CDM-type power spectrum of 
the form 

A 2 m (k) = Q 2 (k/H ) 4 T 2 {k,Q m ,h), (12) 

where we use the CDM transfer function of Bond & Efs- 
tathiou (1984), which gives a reasonable fit for a range of 
CDM models. This provides us with a 3-parameter model 
dependent on the parameter set Q m , h and This is an in- 
teresting parameter set, in particular since CMB data alone 
cannot constrain f2 m without the addition of LSS data. 

If we leave h as a free parameter, the fit diverges off 
to high h and Q and low Q m , preserving the shape, but 
increasing the clustering amplitude. To avoid this we impose 
the HST constraint that h = 0.65 ± 0.12. This effectively 
removes a degree of freedom in the fit. Beyond this we fit 
the PSCz mass points, the abundance of clusters constraint 
and the COBE 4-yr normalisation for a flat universe (Bunn 
& White 1997); 

Q = (1.94 ± 0.2) x 10 _5 fi~ a79 . (13) 

We implicitly assume that the universe is spatially flat, as 
implied by the recent Boomerang (Lange et al 2000) and 
MAXIMA (Hanany et al. 2000, Balbi et al. 2000) results and 
that the spectral index is unity, again as implied by CMB 
results. We also implicitly assume that the contribution to 
the mass density from baryons and neutrinos is negligible. 

Figure 3 shows the \ 2 distribution for the f2 m — Q plane, 
with h — 0.65. The minimum \ 2 values are 

n m = 0.16 ±0.03, Q = (8.4 ±3.8) x 10~ 5 (14) 

where we quote marginal errors. We can also obtain an es- 
timate of the bias parameter for IRAS galaxies: 

6 = 0.84 ±0.28, (15) 

and the spectral shape parameter, 

r = 0.1 ±0.03, (16) 

which is somewhat lower than the usual value of 0.2 from 
the shape of the galaxy clustering spectrum alone. This best 
fit has a x 2 — 2 for 3 degrees of freedom, and so lies within 
the range of acceptable fits. However, while suggestive, this 
analysis makes more assumptions than one would like. These 
assumptions can be dropped by a combined analysis of the 
recent Boomerang and MAXIMA results of the small-angle 
fluctuations in the CMB, which goes beyond our present 
analysis. 

In Figure 2 we plot the PSCz real-space galaxy power 
spectrum and real-space mass power spectrum, along with 



t Our parameter Q is equivalent to the parameter 5jj, the am- 
plitude of clustering at the present Hubble scale, used elsewhere. 
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Figure 3. The x 2 fit to the linear real-space PSCz mass power 
spectrum, the abundance of clusters and the CMB for the two 
parameter space of mass-density parameter, H m , and the ampli- 
tude of mass perturbations, Q. The solid line is the constraint 
from the 4-year COBE analysis of the CMB. 

our best-fit model, a ACDM, normalised to the COBE 4- 
year data. The model fits over two orders of magnitude in 
scale, but nonlinear and nonlocal bias effects lead to a dis- 
agreement on smaller scales, below k = 0.3/iMpc - . Recent 
results from Seljak (2000), Ma & Fry (2000), and Peacock 
& Smith (2000) suggest that this nonlinear regime can best 
be understood in terms of a superposition of randomly po- 
sitioned collapsed clusters, integrated over the cluster mass 
function. Seljak (2000), and Peacock and Smith (2000) go 
further to argue that the nonlinear bias function of galax- 
ies arises due to the statistics of the occupation number of 
galaxies in haloes and show agreement between the ACDM 
model and the small-scale power-law spectrum of galaxies. 
Perhaps the analysis of this regime is not as daunting as it 
once appeared. 



5 CONCLUSIONS 

In this paper we have presented a new analysis of the PSCz 
redshift survey. Using the spherical harmonic analysis of 
HT, BHT and T99 to decompose the survey, taking into 
account the effects of linear redshift-space distortions, non- 
linear "fingers-of-god" effects, limited sky coverage, and the 
radial distribution of galaxies in the survey, we have applied 
a likelihood analysis to a conservative cut of the PSCz sur- 
vey. The catalogue was cut at 0.75Jy and a conservative 
mask used to avoid systematic uncertainties at the low-flux 
end of the catalogue and near the galactic plane. The spheri- 
cal harmonic analysis of the remaining 7042 galaxies resulted 
in 4644 harmonic modes. We used a hierarchical data com- 
pression to reduce this to 2278 for the final analysis. The 
compression was applied using a parameter eigenmode ap- 
proach that compresses information along the line of largest 
degeneracy in parameter space, thus "squeezing" the uncer- 
tainty in this direction. 

We used these methods to estimate the redshift-space 
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distortion parameter, (3 and the amplitude of the real-space 
galaxy power spectrum, parameterised by the amplitude, 
Ao.i, at a wavenumber k = 0.1 /iMpc -1 . Applying the likeli- 
hood analysis to wavenumbers below k = 0.2 /iMpc -1 , where 
linear theory will hold and "fingers-of-god" effects can be 
corrected for, we find 



P = 0.39 ±0.12, 
Ao.i = 0.42 ± 0.02 



(17) 
(18) 



quoting marginalised uncertainties. The distortion param- 
eter is slightly lower than the uncompressed analysis of 
Tadros et al, with an uncertainty reduced by over a fac- 
tor of 2. The consistency of these results with that of the 
earlier analysis of T99 leads us to believe that our analy- 
sis has not been heavily contaminated by nonlinear effects, 
while the conservative cuts in flux and sky coverage avoid 
contamination by flux systematics. These results also are 
in agreement with other, independent determinations of the 
distortion parameter. Comparing the velocity field recon- 
structed from the PSCz with the ENEAR survey, Nusser at 
al (2000) find = 0.5 ± 0.1, while Valentine, Saunders & 
Taylor (2000) find /3 = 0.5 ±0.1 by comparing the recon- 
structed dipole with the observed dipole, the bulk flow of 
galaxies with the MKII survey and the dipole again with 
the SFI catalogue. Finally, using a least squares fit to the 
ratios of the galaxy-velocity to galaxy-galaxy and velocity- 
velocity to galaxy-galaxy power spectra, Hamilton, Tegmark 
& Padmanabhan find (3 = 0.41 ± 0.13 for the PSCz. 

Since in the linear regime the galaxy power spectrum 
is likely to be proportional to the mass power spectrum we 
can combine Ao.i and j3 to find the amplitude of the mass 
perturbations, 



A m (0.1) = (0.16 ±0.04)!^ 



(19) 



Combined with the constraints from the CMB, and HST 
observations of the Hubble parameter, and assuming CDM, 
a flat universe, scale invariant initial mass perturbations 
and negligible contribution to the total mass-density from 
baryons and neutrinos we find Q, m — 0.16 ± 0.03, and a bias 
parameter for IRAS galaxies of b — 0.84 ± 0.28. The mini- 
mum x 2 At nas a value of \ 2 = 2 for 3 degrees of freedom. 
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